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It  cur  initial  pin  of  work,  w  divided  oar  offorta  into  too  aroma: 

■•■•••'•  -v'  -A-'-  - 

1)  Aieljiis  of  thm^fieokfi*  Tech  polar  nephelometer  data;  and 

2)  Analysis  of  the  scattering  data  from  our  nephelometer  and  the  other 
light  acattoring  lnatruaanta  In  relation  to  the  other  aerosol  neaavrenenta 
to  provide  information  that  id  needed  to  model  the  radiative  effects  of  the 
Arctic  base.  Our  final  report,  which  la  structured  in  a  similar  fashion, 
consists  of  a  section  describing  the  AGASP  experiment,  a  section  describ¬ 
ing  the  Geregia  Tech  polar  nephelometer  results,  a  section  describing  oar 
general  analysis  for  modeling  purposes,  and  a  concluding  section. 
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tt»  Arctic  Gm  and  Aerosol  Saapllng  Program  (AGASP)  was  conducts* 


*«*»••  S**  Arctic  in  the  spring  of  19*3  to  study  tits  largo  seals  air 
pollution  phenoasnou  known  as  Aretie  kass.  Tbs  aircraft  nsasuronsnt 
coaponent  of  thlp  progran  utilised^  National  Oceanic  and  Atmospheric 
Adnlnlstratlon  WP-3D  Orion  rsssarch  aircraft  to  study  the  bass  ln-altu  and 
to  ax  tend  into  throe  dimensions  tits  Arctic  bass  asaaursstnts  conducted  at 
tbs  serious  surface  stations  In  tbs  Arctic* 

Tbs  AGASP  prograst  (described  by  Schnell,  1984)  was  a  nultifacstsd 


progran  idtich  Included  tbs  following  aerosol  nsasuranent  objectives:  1. 
Determination  of  aeroOol  else  spOctra,  optical  properties,  chenlcal  conpo- 
sltion,  distribution  and  trajectories}  2.  Measurements  of  both  ln-sttu 
and  surface  observed  radiative  effects  of  the  base.  The  aircraft  aerosol 


Instranentetlon  included  several  aarosol  staples.  a  condensation  nucleus 
counter,  two  PMS  (Particle  Measuring  Syeten)  aerosol  spectroneter  probes 
for  particle  siring,  s  NOAA  Integrating  nephelonster ,  end  the  Georgia  Tech 
polar  nephelonster. 

flight  tracks  of  the  NOAA  VP- 3D  during  the  AGASP  project  are  shown 
in  Figure  2.1  which  Is  taksn  from  Schnell  (1984). 
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Pic*  2.1  Approximate  flight  tracks  of  the  RQAA  WP-3D  during  the  AGASP 
project.  Data  fro*  flights  i  (3-11-83),  2  (3-13-83>, 

3  (3-15-83) ,  and  8  (3-31-83)  ara  discussed  In  this  report 
(fron  Schnell,  1984). 
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Georgia  Tech  polar  nephelometer  was  designed  to  measure  the  total 
aerosol  scattering  and  the  scattering  as  a  function  of  angle,  fro*  which  the 
scattering  phase  function  and  say— etry  parameter  nay  be  deterelned. 

The  basic  concept  of  the  Georgia  Tech  polar  nephelometer  can  be  explained 


with  reference  to  the  device  described  by  Grans  et  al.  (1975).  The  general 
features  of  the  instrument  are  shown  in  Figure  3.1.  A  collimated  laser 


bean  serves  as  the  light  source  and  the  detector  optical  system  defines  a 
narrow  field  of  view  (0.5*  half  angle) .  The  photon  counting  system  amasures 
photomultiplier  pulse  rates  with  the  light  beam  on  and  off,  with  the  differ¬ 
ences  in  the  measured  pulse  rate  being  directly  proportional  to  the  intensity 
of  the  light  scattered  from  the  volume  common  to  the  intersection  of  the 
laser  beam  and  the  detector  field  of  view*  Measurements  are  made  at  diffe¬ 
rent  scattering  angles  by  rotating  the  detector  relative  to  the  direction  of 
propagation  of  the  laser  beam. 

Earlier  Georgia  Tech  measurements  with  polar  nephelometers  have  inclu¬ 
ded  a  ground-based  study  to  determine  the  complex  Index  of  refraction  of 
airborne  soil  particles  (Grams  et  al.,  1974),  a  series  of  airborne  observa¬ 
tions  of  the  optical  properties  of  aerosols  in  the  troposphere  and  lower 
stratosphere  (e.g..  Grams,  1981),  and  laboratory  experiments  on  the  optical 
properties  of  airborne  particles  as  a  function  of  their  size,  shape  and 
refractive  index  (e.g..  Grass  and  Colettl,  1981). 

As  part  of  the  NASA  Aerosol  Climate  Experiment  (ACE) ,  a  new  version  of 
this  Instrument  was  developed  to  be  mounted  and  flown  on  one  of  NASA's  U2 
aircraft.  This  effort  was  Initiated  approximately  5  years  ago  as  part  of 
a  NASA  special  study  which,  in  conjunction  with  satellite  observations  and 
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other  related  aerosol  programs,  was  designed  to  provide  comprehensive  data 
sets  and  analytical  capabilities  for  evaluating  the  impact  of  stratospheric 
aerosol  particles  on  the  Earth's  climate,  for  testing  climate  models  by 
studying  the  properties  and  effects  of  a  large  volcanic  eruption,  and  for 
determining  the  influence  of  aerosols  on  remote  sensing  techniques^ 

the  dew  design  evolved  from  our  efforts  to  satisfy  the  limitations  that 
were  placed  on  the  size,  weight,  and  power  consumption  of  the  instrument 
for  operation  on  the  0-2  aircraft.  The  elimination  of  moving  parts  (such 
as  the  rotating  detector  system  in  the  original  nephelometer)  also  proved 
to  be  highly  desirable  for  aircraft  operations.  We  therefore  designed  a 
system  for  measuring  the  amount  of  light  scattered  by  aerosol  particles 
as  a  function  of  angle  relative  to  the  direction  of  a  collimated  Helium- 
Neon  laser  beam  by  operating  an  image-intensified  scanned-diode  array  in 
the  focal  plane  of  a  wide-angle  "fisheye"  lens.  A  schematic  of  this 
arrangement  is  shown  in  Figure  3.2. 

The  light  scattered  from  the  narrow  collimated  laser  beam  by  atmo¬ 
spheric  molecules  and  aerosol  particles  is  focused  along  a  line  parallel 
to  the  source  beam  in  the  focal  bams  of  the  lens.  The  linear  diode-array 
detector  is  positioned  to  measure  the  light  intensity  along  that  line. 

With  the  180°  wide-angle  lens  focused  at  infinity,  each  point  along  the  array 
corresponds  to  a  discrete  scattering  angle.  The  diode  array  simultaneously 
measures  the  amount  of  light  incident  on  each  of  512  different  positions 
along  the  line.  The  practical  limits  on  the  angular  range  of  the  instru¬ 
ment  are  for  scattering  angl'  from  about  5*  to  175°  from  the  direction 
of  the  source  beam.  In  the  focal  plane  of  the  lens,  this  angular  limit 
transforms  approximately  to  the  middle  400  elements  of  the  diode  array  and 
we  thereby  have  an  angular  resolution  of  approximately  0. 5  * . 


Fig.  3.1.  Schematic  illustration  of  the  Georgia  Tech  laser  polar 
nepheloaeter  (from  Grams •  Daacher ,  and  Wyman,  1975). 
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Fig.  3.2.  Schematic  illustration  of  the  optical  configuration  of  the  Georgia 
Tech  laser  polar  nephelometer  developed  for  operation  on  the  MSA  U-2  aircraft. 
Light  scattered  from  a  laser  beam  by  air  molecules  and  aerosol  particles  is 
focussed  by  a  180°  wide-angle  flsheye  lens  to  a  line  in  the  focal  plane  of  the 
lens  that  is  parallel  to  the  beam  direction.  Light  from  different  angles  is 
focussed  to  different  points  along  that  line.  An  image- intensified  scanned 
linear  diode  array  is  positioned  to  record  the  scattered  light  intensity  at 
each  angle.  In  the  figure,  the  symbol  0j  is  used  to  depict  a  forward  scatter¬ 
ing  angle  and  03  backscatterlng;  02  depicts  light  scattered  normal  to  the  beam 
direction.  See  tert  for  more  details. 


The  signals  recorded  by  all  of  the  elements  of  the  detector  array  are 
periodically  read  by  a  microprocessor  based  data  management  system.  This 
system  controls  the  operation  of  the  Instrument  and  records  data  using  a 
Digital  Equipment  Corporation  LSI-11/23  microprocessor  with  64  Kbytes  of 
memory  with  a  disk  operating  system  based  on  a  19  Mbyte  Winchester  disk 
drive.  The  disk  recording  system  is  used  to  store  the  instrument  control 
and  data  analysis  programs  and,  also,  to  record  the  large  amounts  of  data 
read  from  the  detector  array. 

The  laser  beam  used  in  the  nephelometer  is  linearly  polarized.  A  disk 
in  the  path  of  the  laser  beam  can  be  rotated  to  3  different  positions  by  a 
small  motor.  One  position  blocks  the  beam  so  that  a  background  signal  without 
the  presence  of  the  source  beam  can  be  established.  The  second  position  is  a 
circular  hole  in  the  disk  which  allows  the  linearly  polarized  laser  beam  to 
illuminate  the  volume  being  sampled;  the  location  of  the  detector  system  is 
designed  to  allow  us  to  observe  the  polar  dlstrlbuiton  of  light  scattered  in 
the  plane  that  is  perpendicular  to  the  direction  of  the  electric  vector  of 
the  source  beam.  A  half-wave  plate  has  been  mounted  at  the  third  position; 
this  optical  element  rotates  the  original  plane  of  polarization  of  the  laser 
beam  by  90  degrees  so  that  the  detector  system  measures  the  polar  diagram 
for  the  parallel  scattering  plane.  Each  nephelometer  scan  consists  of  the 
signals  recorded  by  each  of  the  512  elements  of  the  detector  array  for  the 
3  different  positions  of  the  disk.  The  sample  time  required  to  complete  each 
3-part  scan  can  vary  from  just  a  few  seconds  for  use  in  conditions  of  heavy 
aerosol  loading  to  longer  time  intervals  for  observations  in  high  visibility 
conditions.  To  indicate  the  overall  sensitivity  of  the  system,  we  point  out 
that  integration  times  of  approximately  one  minute  are  required  to  obtain 
observations  of  molecular  scattering  (without  aerosol  particles  present) 
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when  operating  the  Instrument  at  U-2  sampling  altitudes  with  ambient  pressures 
of  about  100  millibars. 

The  construction  of  the  nephelometer  was  completed  and  the  device  was 
first  flown  in  test  flights  on  the  NASA  U-2  during  October  1981  and  April 
1982.  Our  first  data  flights  were  made  during  July  1983.  Additional  strato¬ 
spheric  data  were  obtained  during  U-2  and  ER-2  flights  made  in  December  1982 
and  July  1983.  With  respect  ot  the  current  study,  the  nephelometer  was  also 
flown  as  part  of  the  aerosol  instrument  package  on  the  NOAA  P-3  aircraft  as 
part  of  the  Arctic  Gas  and  Aerosol  Sampling  Program  (AGASP)  study  of  Arctic 
haze  layers  (Schnell,  1984). 

The  potential  for  generating  large  amounts  of  data  during  field  measure¬ 
ments  that  involve  several  6-8  hour  flights  is  clearly  present  with  the 
instrument.  For  example,  our  U-2  flights  provided  us  with  almost  850  three 
part  scans  during  the  July  1982  flights,  approximately  300  scans  during 
the  December  1982  flights  and  several  hundred  additional  scans  during  the 
July  1983  flights.  The  Arctic  haze  observations  during  the  AGASP  flights 
in  Spring  1983  provided  almost  3000  more  3-part  scans.  A  significant  part 
of  our  data  analysis  effort  has  involved  the  development  of  procedures  to 
present  such  large  amounts  of  raw  data  in  scientifically  useful  forms. 

Our  initial  approach  was  that  of  preparing  graphs  or  tables  of  scatter¬ 
ing  cross  sections  and  scattering  asymmetry  parameters  obtained  by  applying 
least-squares  curve-fitting  techniques  to  establish  "best-fit"  parameters 
for  Henyey-Greenstein  phase  functions  as  described  by  Grams  (1981).  In 
support  of  qualitative  "quick-look"  applications,  we  also  carried  out  some 
analytical  studies  showing  that  the  amount  of  light  detected  by  the  polar 
nephelometer  at  scattering  angles  near  45°  is  approximately  proportional 
(to  within  a  factor  of  2)  to  the  total  scattering  cross  section  of  the 
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aerosol  particles*  An  example  of  a  plot  of  45°  scattering  cross  section 
versus  time  as  recorded  by  the  laser  nephelometer  on  a  U-2  flight  is  shown 
in  Figure  3*3*  Data  of  this  type  were  printed  out  in  real  time  on  the 
NQAA  P-3  aircraft  during  the  Spring  1983  AGASP  flights  to  help  establish 
the  times  that  the  aircraft  was  within  a  layer  of  Arctic  haze. 

During  the  past  year,  we  Initiated  a  comprehensive  program  for  cali¬ 
brating  the  nephelometer.  Although  the  recording  of  data  with  the  instru¬ 
ment  could  be  (and  was)  accomplished  without  an  absolute  calibration  of  the 
instrument,  we  felt  that  our  results  could  not  be  published  until  such 
calibrations  had  been  completed.  In  some  of  our  field  observations,  such 
as  the  AGASP  Arctic  haze  flights,  the  real-time  45°  output  of  the  instru¬ 
ment,  as  illustrated  in  Figure  3.3,  was  used  by  the  scientists  aboard  the 
aircraft  as  a  means  of  determining  when  the  aircraft  was  within  a  layer  of 
haze  particles.  However,  one  of  our  major  goals  has  been  to  obtain  a  more 
quantitative  treatment  of  the  data. 

Perhaps  the  most  difficult  factor  in  relating  the  photon  counts 
recorded  by  each  channel  of  the  diode  array  has  been  the  extreme  diffi¬ 
culty  in  establishing  the  exact  proportion  of  the  observed  counts  due  to 
molecular  scattering  by  the  atmospheric  volume  being  sampled.  To  illus¬ 
trate  the  nature  of  the  problem,  the  reader  is  referred  to  Figure  3.4 
which  is  a  published  example  of  a  polar  nephelometer  observation  made  by 
an  earlier  version  of  the  Georgia  Tech  laser  polar  nephelometer  system 
(Grams,  1981) .  While  the  figure  refers  to  observations  of  light  scattering 
by  aerosol  particles  in  the  stratosphere,  it  is  representative  of  the  type 
of  observations  provided  by  our  airborne  nephelometer  system. 

The  curves  shown  in  the  figure  illustrate  several  important  points. 

The  first  and  most  Important  point  to  be  made  is  the  fact  that,  except  for 
the  extreme  forward  directions,  molecular  scattering  often  dominates  the 


9 


f 


1 

*S* 

742 

5 

S30 

* 

MS 

s 

J99 

• 

•33 

HtlMtlMMHIHUimn 

rm 

•*•*••*** ««**«9«*t *««*■« 

.  • 

45* 

949*994491194 

4 

394 

4449*4***** 

to 

377 

9*49949**9* 

OOl 21 |04  12/17/03 

tt 

224 

**99**9* 

12 

317 

13 

30* 

14 

32S 

19 

24* 

*99*99* 

1* 

314 

17 

390 

IS 

313 

********* 

10 

290 

30 

293 

OOl 35 <37  12/17/02 

31 

222 

99*99999 

mm 

IS* 

9***9 

33 

20* 

94*9*9 

34 

241 

******* 

39 

1  *7 

mu 

3* 

19* 

imn 

37 

13S 

*99* 

39 

137 

*9* 

39 

13S 

•**• 

30 

13S 

•4*9 

00.43.40  12/17/02 

21 

1*9 

***• 

33 

93 

II 

»  + 

99 

•  •  • 

:« 

90 

«*• 

39 

119 

IK 

24 

111 

*44 

27 

112 

•99* 

31 

111 

*4* 

39 

109 

•*• 

*0 

123 

IK 

*>•321 17  12/17/32  " 

41 

107 

**• 

*2 

120 

•4* 

4 Z 

04 

•9 

44 

lOl 

*99 

49 

77 

*• 

44 

107 

IK 

47 

99 

49* 

4S 

M 

*4 

49 

US 

*•• 

30 

9* 

•9* 

01*04*04  12/ 17/02 

21 

91 

•  t 

23 

lOl 

*99 

22 

*3 

*9 

14 

102 

III 

32 

*0 

4 

34 

*4 

II 

27 

•• 

II 

SS 

31 

4 

24 

•n 

•  4 

40 

9* 

4 

Fig.  3.3.  Examples  of  low  resolution 
output  of  45*  scattering  cross 
section  versus  tins  as  recorded  by 
the  laser  naphalometer  system  as  the 
NASA  U-2  aircraft  made  a  vertical 
ascent  during  one  of  the  U-2  flights 
for  NASA's  Aerosol  Climate  Experiment 
during  December  1983.  As  pointed  out 
In  the  text*  this  parameter  is 
approximately  proportional  to  the 
total  scattering  cross  section. 
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Fig.  3.4.  Exanple  of  polar  nephelonetar  observations  for  Che 
stratospheric  aerosol  layer  over  Poker  Flat.  Alaska.  The  data 
were  obtained  on  July  17,  1979.  The  squares  are  data  points  for 
scattering  In  the  plane  perpendicular  to  the  electric  vector  of 
Incident  polarised  light;  the  circles  are  data  points  for  the 
parallel  plane.  Dotted  lines  are  relative  contributions  to  the 
phase  functions  due  to  Molecular  scattering;  solid  lines 
represent  the  phase  function  for  both  aerosol  and  nolecular 
scattering  using  the  Banyey-Greeosteln  phase  function  to  specify 
the  aerosol  contributions.  Best-fit  Henyey-Greenstein  parameters 
are  an  asyunetry  parameter  g  -  0.32  end  a  ratio  of  aerosol  to 
nolecular  scattering  coefficient  f  -  0.65. 


observed  light-scattering  pattern.  Thus*  relatively  snail  errors  in  remov¬ 
ing  the  aolecvlar  scattering  contributions  can  seriously  affect  the  inferen¬ 


ces  that  one  uakes  about  the  aaount  and  angular  distribution  of  light 
scattered  by  airborne  particles.  The  second  point  to  be  made  is  the  fact 
that  the  observed  scattering  patterns  can  usually  be  fit  vith  the  so-called 
Henyey-Greensteln  phase  function  (Grans*  1981) 

PHG(0)  "  (1  +  g*(-  2g*cos  6)-*/^ 

in  which  g  la  the  scattering  asyneetry  parameter  and  9  is  the  scattering 
angle.  Figure  3.5  shows  the  fora  of  these  curves  (dotted  lines)  for  3 
cases  ranging  froei  g  *  0.5  (botton  curve)  to  g  ■  0.7  (top  curve).  While 
the  procedure  described  by  Grans  (1981)  for  obtaining  "best-fit"  Henyey- 
Greensteln  parameters  given  reasonable  results  for  describing  the  scatter¬ 
ing  pattern  for  the  sun  of  the  aerosol  and  molecular  components  for  cases 
in  which  the  aerosols  are  relatively  small  contributors  to  the  scattered 
light*  cases  involving  heavy  aerosol  loading  can  cause  complications  for 
the  analysis  procedures  and  we  found  it  necessary  to  search  for  a  better 
approach . 

At  the  present  time,  we  believe  that  a  modified  version  of  the  tech¬ 
nique  described  by  Grans  (1981)  will  provide  significantly  better  aerosol 
scattering  cross  sections.  This  technique  takes  advantage  of  the  wide 
range  of  angles  measured  with  the  nephelometer  device  to  provide  a  parameter 
that*  in  principle*  completely  eliminates  contributions  due  to  molecules. 
This  procedure  Involves  calculating  for  each  forward  angle  the  difference 
between  the  observation  at  that  angle  and  the  corresponding  observation  at 
the  complementary  angle.  Thus*  since  the  molecular  scattering  phase 
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SCATTERING  ANGLE 


Fig.  3.5.  Dotted  lines  srs  the  Hsnysy-Grssnstsin  phase  functions  for 
g  -  0.5  (bottoa  line)  to  g  -  0.7  (top  line).  Solid  lines 
srs  the  results  of  subtracting  phase  function  values  at 
coapleaentary  angles  froa  values  of  the  Hsnysy-Grssnstsin 
phase  functions  in  the  forward  directions. 


function  Is  symmetric  about  90*,  the  subtraction  will  remove  all  aolecular 
contrlbutlona  to  the  nodlfiad  scattering  parameter.  The  differences  compu¬ 
ted  at  each  angle  will  than  be  directly  proportional  to  aerosol  concentra¬ 
tions  and  independent  of  the  local  aolecular  number  density.  The  solid 
lines  in  Fig.  3.5  show  the  results  of  computing  such  differences  for  the 
Henyay-Greensteln  phase  functions.  Note  that  for  angles  greater  than  about 
40*,  the  curves  are  almost  independent  of  the  choice  of  the  asymmetry 
parameter  value.  Our  preliminary  analysis,  used  to  generate  the  scatter¬ 
ing  cross  sections  presented  in  this  report,  applies  our  least-squares 
curve-fitting  programs  to  the  problem  of  fitting  the  observations  to  the 
modified  Henyey-Greensteln  functions  for  angles  between  40  and  90*, 
assuming  g  ■  0.6  (which  is  typical  of  the  "best-fit"  asymmetry  parameters 
obtained  from  our  AGASP  observations).  This  approach  is  a  relatively 
recent  development  and  further  study  is  warranted.  However,  our  prelimi¬ 
nary  analysis  involving  the  use  of  Henyey-Greensteln  phase  functions 
Indicates  that  the  absolute  values  of  the  scattering  cross-sections 
obtained  by  this  approach  should  be  accurate  to  within  about  20Z. 

Figures  3.6  and  3.8  show  the  results  of  our  analysis  for  two  days  of 
AGASP  data  -  March  21  and  March  31,  1983.  Figures  3.7  and  3.9  show  the 
corresponding  analyses  for  condensation  nuclei,  oxone,  aerosol  extinction 
(integrating  nephelometer) ,  wind,  temperature  and  dew  point  data  obtained 
by  other  investigators  on  the  same  flights  (Raatz  et  al.,  1985).  In  subse¬ 
quent  discussions,  we  point  out  that  the  integrating  nephelometer  results 
are  almost  an  order  of  magnitude  higher  than  those  calculated  from  the 
size  distributions  obtained  with  the  particle  sizing  Instruments  on  the 
AGASP  flights.  Our  first  comment  on  comparing  our  results  with  the  other 
instruments  is  that  the  scattering  cross  sections  obtained  from  the  polar 
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nephelometer  data  are  closer  to  those  measured  with  integrating  nephelometer. 
In  most  cases*  the  polar  nephelometer  values  are  somewhat  smaller  than 
those  measured  with  the  integrating  nephelometer.  However,  in  situations 
in  which  a  high  amount  of  variability  in  aerosol  concentration  is  evident, 
the  polar  nephelometer  appears  to  respond  faster,  increasing  to  relatively 
higher  values  than  those  recorded  by  the  integrating  nephelometer.  The 
two  different  types  of  nephelometers  do  seem,  therefore,  to  be  in  general 
agreement ,  although  detailed  comparisons  often  do  tend  to  lead  to  some 
interesting  differences.  For  example,  at  about  1810Z  on  Marfeh  21,  Figs. 

3.7  and  3.8  show  a  case  in  which  the  normal  situation,  as  described  above, 
is  reversed;  namely,  the  integrating  nephelometer  picked  up  a  significant 
increase  in  aerosol  extinction  which  was  just  barely  detected  by  our 
instrument. 
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.  GENERAL  ANALYSIS  OF  SIZE  AND  SCATTERING  DATA 


This  analysis  was  primarily  concerned  with  determination  o£  best 


estimates  of  particle  size  for  use  in  models  of  the  radiative  properties 


of  the  aerosol  as  well  as  for  the  interpretation  of  the  scattering  data. 


Among  the  areas  of  Interest  were  the  testing  of  the  hypotheses  incorpora¬ 


ted  in  the  Patterson  et  al.  (1982)  radiative  models  and  the  estimation  of 


the  radiative  properties  at  infrared  wavelengths.  We  have  also  identified 


some  significant  inconsistencies  between  two  of  the'  AGASP  data  sets. 


Although  not  an  original  goal,  we  have  attempted  to  reconcile  these  two 


data  sets,  those  of  the  integrating  nephelometer  and  the  PMS  optical 


particle  counters,  and  to  understand  the  reasons  for  the  inconsistencies. 


Size  Distributions  -  General  Comparisons 


The  AGASP  flights,  other  aircraft  programs,  and  ground-based  measure¬ 


ment  programs  in  1983  provided  many  different  measurements  of  the  aerosol 


size  distribution  under  a  variety  of  haze  conditions.  Consideration  of  a 


representative  set  of  these  measurements  together  with  some  earlier  size 


distribution  data  can  provide  useful  insight  Into  the  appropriate  mode 


parameters  for  use  in  radiative  modeling  of  the  Arctic  haze. 


A  comparison  of  the  data  is  shown  in  Figure  4.1,  which  is  a  rather 


complex  summary  of  the  different  size  distributions.  Each  of  the  curves 


on  this  figure  represent  published  data  from  measurements  made  during  a 


spring  aerosol  maximum  -  either  1983  or  a  comparable  earlier  period  (the 


Belntzenberg  data).  The  published  data  have  all  been  converted  to 


dV/d  log  D  distributions  for  consistency  of  presentation. 


6/**err<  (4/m) 


Size  distribution  measurements  of  Arctic  haze  aerosols  from 
a  variety  of  sources  plotted  as  dv/d  log  d  distributions. 


Two  of  the  curves  in  the  figure  represent  size  distributions  measured  during 

AGASP  flights (3-1 1-83,  - ;  3-31-83, - )  [Schnell  and  Raatz,  1984; 

Raatz  and  Schnell,  1984].  The  other  curves  show  a  size  distribution  mea¬ 
sured  over  the  Arctic  Ocean  by  the  University  of  Washington  group  (••••) 

[Radke  et  al.,  1984];  two  1983  size  distributions  measured  at  surface  sites 
in  Alaska,  one  at  Barrow  (o-o-o;  Bodhaine,  1984)  and  one  in  central  Alaska 
(*-*-*;  Shaw,  1984);  a  size  distribution  measured  at  Spitzbergen  by  Heint- 
zenberg  (1980)  (-+-+-+) ;  and  a  size  distribution  measured  in  the  Canadian 
Arctic  ( - )  by  Huff  and  Trivett. 

We  have  also  shown  for  comparison  portions  of  the  two  modes  used  by 
Patterson  et  al.  (1982)  for  their  model  of  the  Barrow  spring  maximum. 

These  models  of  Patterson  et  al  are  shown  displaced  from  the  other  modes 
for  clarity. 

Figure  4.1  shows  quite  clearly  the  overall  modal  structure  of  the 
Arctic  hazfe  aerosol.  Each  of  the  distributions  which  covers  the  submicron 
particle  size  range  shows  a  submicron  particle  mode  volume  maximum  between 
v  0.2  and  0.3  ym  diameter.  Each  of  the  distributions  which  covers  the 
giant  particle  size  range  shows  a  peak  in  the  2  to  3  ym  size  range. 

The  0.3  ym  mode,  which  consists  of  sulfuric  acid  droplets  as  well  as 
other  secondary  and  combustion  aerosols,  is  quite  similar  in  mode 
characteristics  to  the  accumulation  mode  measured  in  mid-latitude  regions 
as  discussed  by  Patterson  et  al.  (1982)  and  to  measurements  made  using 
microscopic  analysis  techniques  (Biggs,  1980).  The  consistency  in  the 
different  data  sets  for  this  accumulation  mode  indicates  that  we  have  a 
reasonable  knowledge  of  the  mode  parameters  of  this  major  peak.  Mie 
calculations  indicate  that  the  differences  in  the  modes  will  not  be  signi¬ 
ficant  when  the  calculations  are  normalized  to  measured  scattering  or 
extinction. 
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Thare  is  somewhat  greater  variability  in  observed  mode  parameters  of 
the  supermicron  giant  particle  mode.,  but  all  of  the  data  sets  are  consis¬ 
tent  in  indicating  a  mode  volume  maximum  between  1.8  and  4  um  diameter. 

For  these  particles,  the  Patterson  study  inferred  a  mean  radius  somewhat 
larger  than  was  generally  observed  during  the  AGASP  flights,  presumably 
due  to  longer  transit  times  for  these  Arctic  aerosols.  Microscopic 
analysis  (Rahn.  1978;  Bailey  et  al.,  1984)  shows  that  giant  particles 
consist  in  large  part  of  crustal  particles;  ionic  composition  data  of  Radlte 
(1984)  suggest  that  there  can  be  a  significant  sulfate  and  sea  salt  compo¬ 
nent  for  these  larger  particles  as  well.  The  exact  composition  of  this 
giant  particle  mode  will  vary  with  location  and  time. 

There  is  an  additional  feature  in  three  of  the  data  sets,  those  of 
Shaw,  Schnell,  and  Radke,  which  is  a  mode  or  a  shoulder  near  0.7  pm 
diameter.  The  reason  for  this  feature  is  not  known,  although  the  presence 
of  the  feature  in  three  separate  data  sets  suggests  that  it  is  not  an 
artifact.  In  fact,  Hoppel  et  al.  (1985)  report  the  presence  of  a  similar 
mode  in  marine  boundary  layer  data  and  discusses  possible  production 
mechanisms  for  such  a  mode. 

While  there  is  variation  in  the  magnitudes  of  the  modes  and  in  the 
presence  or  absence  of  a  0.7  im  diameter  particle  mode,  we  can  say  that, 
in  general,  the  broad  features  of  the  size  distribution  for  the  spring 
aerosol  appear  to  be  surprisingly  well  defined  between  0.1  and  10  ym 
diameter.  In  addition,  the  similarity  of  these  features  to  mid-latitude 
features  confirms  that  the  overall  processes  determining  the  size  distri¬ 
bution  can  be  understood  in  terms  of  similar  mid-latitude  processes,  con¬ 
firming  the  inferences  of  Patterson  et  al. 


Size  Distributions  -  AGASP  Data 

Because  of  the  ability  to  make  detailed  comparisons  with  the  optical 
and  the  chemical  data,  including  the  integrating  and  polar  nephelometer 
data,  we  have  emphasized  interpretation  of  the  size  data  obtained  during 
the  AGASP  program. 

The  real-time  size  distribution  measurements  aboard  the  WP-3D  during 
the  AGASP  flights  were  made  with  two  externally  mounted  PMS  aerosol  probes, 
an  ASASP-100X  for  smaller  particles  (those  with  d  <  3  ym)  and  an  FSSP  for 
larger  "giant"  particles  (those  with  diameters  between  3  and  10  ym) . 
Analysis  has  shown  that  both  probes  are  needed  to  measure  particle  size 
for  radiatively  important  particles. 

Six  representative  measurements  of  the  aerorol  size  distributions  are 
shown  in  Figs.  4.2a,  4.2b,  and  4.2c.  The  figures  include  ASASP-100X 
data  (•)  and  FSSP  data  (o) ;  the  presentation  is  in  the  form  of  number  con¬ 
centration  as  well  as  mass  concentration  versus  size.  The  plots  in  Figs. 
4.2a  and  4.2b  represent  data  in  the  Alaskan  Arctic;  the  plots  in  Fig.  4.2c 
represent  data  in  the  European  Arctic. 

It  is  apparent  that  in  these,  as  in  the  other  data  sets,  the  aerosol 
exhibits  a  blmodal  structure.  The  mode  with  a  mass  peak  near  2  ym  is  well 
defined.  Because  of  some  problems  with  the  ASASP  probe,  the  first  three 
channels  of  data  were  not  used  in  the  analysis.  Consequently,  the  data 
extend  only  to  0.3  ym  diameter  and  are  not  sufficient  to  determine  mode 
parameters  for  the  submicron  mode.  On  the  basis  of  aerosol  microphysics 
as  well  as  the  other  Arctic  measurements  discussed  above,  we  can  infer 
that  the  mass  peak  is  near  0.3  ym  diameter.  This  assumption  is  not  incon¬ 
sistent  with  the  PMS  optical  particle  counter  data. 
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Previous  work  has  shown  that  aerosol  size  distributions  can  generally 


be  represented  by  a  log-normal  distribution  of  the  form 


dW  _  Sq 
d  log  r  /2ir  log  cf 


^Clog  r  -  log  ro)2j 
2log  2o 


in  which  NQ  is  the  total  particle  number.  r0  is  the  number  mean  radius,  and 
or  is  the  parameter,  the  standard  deviation,  which  describes  the  width  of 
the  distribution.  For  the  log-normal  distribution,  the  volume  distribution 
has  the  same  form  as  the  number  distribution,  and  the  volume  mean  radius  is 
related  to  the  number  mean  radius  by  the  equation 


In  rv  ■  In  tN  +  3  In2 a 


The  mass  distribution  may  be  determined  by  multiplying  the  volume  distri¬ 
bution  by  the  appropriate  particle  density.  The  radius  terms  may  be 
replaced  with  diameter  terms  in  either  of  the  equations. 

We  have  made  approximate  log-normal  fits  to  the  giant  particle  mode 
centered  near  3  pm  on  the  volume  distributions.  Although  there  were 
minor  differences  in  mode  parameters  in  each  case,  an  adequate  fit  for  all  of 
the  distributions  was  obtained  with  d^  value  of  1  ym  and  a  a  of  1.75,  with 
a  corresponding  dy  value  of  2.6  ym.  We  see  that  both  the  European  Arctic 
and  the  Alaskan  Arctic  data  show  a  striking  similarity  in  mode  parameters. 
This  suggests  that  these  giant  particles  are  representative  of  a  well  aged 
aerosol  that  has  been  transported  a  long  distance. 

The  data  suggest  a  good  deal  of  variability  in  the  submicron  particle 
mode.  The  data  as  presented  in  Figs.  4.2a,  b,  c  suggest  that  the  giant 
particles  are  relatively  more  important  in  producing  the  radiative  effects 
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of  the  aerosol  that  w*  would  expect  on  the  heels  of  other  dete  sets.  This 
Is  not  the  case,  however,  as  we  will  see  in  the  next  section,  since  a  wore 
detailed  analysis  of  the  PMS  data  set  Indicated  some  problems  in  the  data. 


particularly  with  regard  to  the  submicron  mode  concentration  and  sizing. 
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The  validity  of  an  Individual  data  set  may  be  tested  by  means  of  inter- 
comparisons  among  simultaneous  data  sets.  The  large  number  of  simultaneous 
measurements  of  aerosol  parameters  measured  in  AGASP  allows  us  to  make 
such  detailed  intercomparisons  among  the  measurement  sets.  One  such  inter- 
comparison  involves  the  PMS  data  and  the  simultaneous  integrating  nephelo- 
meter  data.  In  this  comparison  the  PMS  size  distributions  presented  by 
Raatz  and  Schnell  were  used  to  calculate  the  scattering  expected  for  the 
aerosol  for  a  wavelength  of  0.55  wm,  the  approximate  wavelength  of  the 
integrating  nephelometer.  Separate  calculations  were  made  for  the  sub- 
micron  and  the  giant  particle  modes;  nominal  optical  properties  (Patterson 
et  al.,  1982)  were  used  for  both  the  giant  and  the  submicron  particle 
mode  calculations.  As  a  first  cut,  because  of  the  uncertainties  of  fitting 
a  log-normal  distribution  to  the  submicron  PMS  data,  this  submicron  data 
was  approximated  by  a  Junge  distribution  extending  from  a  diameter  of 
0.2  pm  to  1  ym.  The  results  of  this  calculation  for  scattering  are  shown 
in  Table  4.1;  simultaneous  integrating  nephelometer  data  are  also  shown 
in  Table  4.1. 

It  is  apparent  from  the  data  in  Table  4.1  that  there  is  a  serious 
discrepancy  between  the  optical  scattering  inferred  from  the  AGASP  PMS 
data  and  the  nephelometer  scattering  data,  with  the  calculated  scattering 
less  than  the  observed  scattering  by  an  order  of  magnitude  or  more.  The  polar 
nephelometer  data,  although  somewhat  lower  than  the  integrating  nephelometer 
data,  show  the  same  discrepancy. 


■>53W 


■yi/f. 


TABLE  4.1 

PMS  DISTRIBUTION  SCATTERING  CALCULATIONS 


DATE 


HAZE 

CONDITION 


SCATTERING  CALCULATIONS  FOR 
X  ■  .55  in*  (a”1) 


Submlcron 


Giant 


Total 


H  «  Haze.  C  *  Clear 
Schnell  and  Raatz  (1984) 
Raatz  and  Schnell  (1984) 
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TABLE  4.2 


DATE 


MODIFIED  SCATTERING  CALCULATIONS 


HAZE 

CONDITIONS 


SCATTERING  CALCULATIONS  FOR 
X  -  .55  pa  (a"1) 


Submicron 


Giant 


Total 


H  -  Haze.  C  -  Clear 
1Schnell  and  Raatz  (1984) 
2Raatz  and  Schnell  (1984) 
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TABLE  4.3 


10.6  vm  EXTINCTION  DATA 

DATE  RAZE  CONDITION  CALCULATED  EXTINCTION  (nf1) 


Submicron 

Giant 

Total 

3-11-83 

Haze 

2.4  x 

io"7 

2.1  x 

10~7 

4.5  x  10” 

3-11-83 

Clear 

.5  x 

io"7 

.4  x 

K 

1 

o 

r4 

0.9  x  10“ 

3-13-83 

Haze 

5.4  x 

io"7 

.8  x 

io”7 

6.2  x  10“ 

3-15-83 

Haze 

2.1  x 

io"7 

4.0  x 

io"7 

6.1  x  10" 

3-31-83 

Haze 

3  x 

io'7 

6.9  x 

10"7 

9.9  x  10" 

3-31-83 

Clear 

.9  x 

io"7 

.5  x 

10’7 

1.4  x  10" 
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The  Integrating  nepheloaeter  date  may  also  be  comapred  with  the 
direct  radiative  mesurements  of  Valero  et  al.  (1983)  for  March  31* 1983. 
Valero  aeasured  upward  and  downward  flux  at  several  wavelengths.  Ackerman 
and  Valero  (1984)  used  this  flux  data  and  assumed  nominal  aerosol  proper¬ 
ties  to  estimate  an  optical  depth  of  \  0.3  for  the  haze  layers  at  mid- 
visible  wavelengths.  By  comparison,  an  approximate  Integration  of  the 
integrating  nephelometer  data  with  altitude  suggests  an  optical  depth  due 
to  scattering  of  about  .13,  roughly  half  of  the  total  optical  depth 
Inferred  by  Ackerman  and  Valero.  Although  the  numbers  are  not  entirely 
consistent,  the  comparison  does  indicate  that  the  nephelometer  is  measuring 
scattering  that  is  not  grossly  in  error. 

A  further  comparison  of  observed  volume  concentrations  of  the  submicron 
particle  mode  with  mass  concentrations  of  carbon  and  sulfur  (which  are 
primarily  associated  wtih  the  submicron  mode)  presented  by  Hanson  and 
Rosen,  1984  for  the  March  11  and  March  13  flights  indicates  that  the  PMS 
data  significantly  underestimates  the  amount  of  material  present  in  the 
small  particle  mode,  possibly  by  a  factor  of  4  or  more.  These  compari¬ 
sons  suggest  very  clearly  that  the  PMS  data  represents  a  significant 
underestimation  of  the  actual  aerosol  concentrations. 

The  lack  of  small  particle  data  and  the  possibility  of  response 
differences  in  the  PMS  probes  between  the  calibration  aerosol  and  the 
measured  aerosol  (see  e.g.,  Pinnick  et  al.,  1983  for  discussion)  in  the 
0.3  -  0.5  urn  size  range  suggest  that  the  discrepancy  may  be  primarily  in 
the  smaller  size  particles.  We  have  tested  this  possibility  by  assuming 
that  the  actual  submicron  size  distribution  is  an  average  size  distribu¬ 
tion  that  is  determined  from  the  curves  on  Fig.  4.1.  We  have  used  the 
Shaw  data  from  the  figure  to  fit  a  log-normal  distribution  with  a  dg  of 


0.17  ua  a  dy  of  .28  pa  and  a  a  of  1.46.  The  appropriate  particle  number 
for  the  dlatrlbutlon  Is  assuaed  to  be  given  by  the  total  particle  number 
aeaaured  with  the  condensation  nucleus  counter  aboard  the  VP- 3D  aircraft 
for  the  tlae  period  In  question. 

The  results  of  these  new  calculations  are  shown  In  Table  4.2.  It  Is 
apparent  that  the  discrepancy  between  calculated  and  measured  scattering 
has  been  significantly  reduced.  This  Indicates  that  the  assuaption  that 
the  PMS  probes  underestimate  the  saall  particle  concentration  is  a  valid 
one.  We  Infer  that  a  size  distribution  determined  by  combining  an  average 
submicron  size  distribution  that  is  matched  to  the  measured  particle 
number  with  the  aeasured  giant  particle  size  distribution  Is  adequate  as 
a  first  approximation  to  the  actual  size  distribution. 

An  example  of  such  a  modified  size  distribution  in  shorn  in  Figure 
4.3.  We  can  see  In  the  figure  that  the  PMS  optical  particle  counter  data 
appear  to  underestimate  our  Inferred  particle  number  in  the  0.3  -  0.5  ua 
size  range.  This  aay  be  explained  with  reference  to  data  of  Pinnlck  and 
Auveraann  (1979)  who  have  shown  that  the  ASASP  probes  will  underestimate 
the  particle  size  in  this  size  range  if  the  particle  is  highly  absorbing 
relative  to  the  calibration  aerosol.  The  Arctic  aerosol  is  known  to  be 
absorbing  and  the  real  component  of  the  refractive  index  is  less  than  that 
of  the  polystyrene  latex  spheres.  We  would  expect  then  that  the  measured 
particles  are  larger  than  indicated ,  which  would  have  the  effect  of  increas' 
ing  the  aaxlaua  size,  reducing  the  logarithmic  width  of  the  size  range, 
and  causing  the  dV/d  log  d  plot  to  be  higher  than  indicated.  This  would 
result  In  a  greater  consistency  between  the  Inferred  and  the  measured 
distributions. 
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Inferred  giant  and  submicron  scattering  for  the  AGASP  data 
(•)  and  for  the  Radke  et  al.  (1984)  haze  data.  The  figure 
confine  the  lnportance  of  the  submicron  mode. 


Inferences  from  Size  and  Scattering  Data 


The  purpose  of  this  section  Is  to  use  our  size  data,  as  modified  to 
conform  with  the  scattering  data,  to  give  an  estimate  of  the  relative 
Importance  of  submicron  and  giant  particles  in  producing  visible  wavelength 
radiative  effects  and  to  determine  the  extinction  of  the  aerosol  at  10.6  am 
as  an  estimation  of  the  magnitude  of  the  infrared  effects. 

We  have  used  our  modified  size  distribution  to  calculate  scattering 
for  both  submicron  and  giant  particles.  The  results  for  each  of  the  AGASP 
samples,  and  for  the  Radke  data,  are  shown  in  Fig.  4.4a.  The  modified 
AGASP  data  is  consistent  with  the  Radke  data  in  confirming  that  the  giant 
particle  mode  is  of  little  importance  in  producing  visible  wavelength 
effects  in  haze  situations. 

The  mid-infrared  effects  are  expected  to  be  much  more  heavily 
influenced  by  the  giant  parltcle  mode  and  this  is  confirmed  by  the  data  of 
Table  4.3  which  shows  the  extinction  of  10.6  ym  due  to  both  submicron  and 
giant  particle  modes  for  each  of  the  data  sets,  as  well  as  the  total  extinc¬ 
tion.  These  data  were  calculated  using  the  urban  aerosol  refractive  index 
data  of  Shettle  and  Fell  (1979)  for  the  submicron  particles;  the  giant 
particle  refractive  index  was  assumed  to  be  1.7  -  .31,  a  value  appropriate 
to  soil  particles  as  discussed  by  Patterson  (1981). 

The  calculated  values  for  the  10.6  ym  extinction  in  haze  layers  range 
from  .4  x  10-7m_1  to  0.9  x  10~7m_1;  possible  variations  in  the  submicron 
particle  concentration  will  have  an  effect  on  the  calculated  extinction, 
but  the  effect  is  less  significant  than  at  visible  wavelengths.  Those 
10.6  ym  extinction  values  are  somewhat  more  than  an  order  of  magnitude 
less  than  the  visible  wavelength  values.  The  Infrared  effects  are  of 
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potential  importance,  however,  because  the  effects  can  persist  throughout 
the  polar  winter  night  while  the  visible  wavelength  effects  are  of  impor¬ 
tance  only  in  a  relatively  short  period  of  time  in  the  spring.  The  possible 
importance  of  these  infrared  effects  will  have  to  be  determined  by  appro¬ 
priate  modeling  efforts. 

Because  of  the  importance  of  the  absorption  in  determining  the  infrared 
effects,  knowledge  of  the  imaginary  component  of  the  refractive  index  is 
crucial  in  determining  these  effects.  There  is,  however,  considerable 
uncertainty  in  the  appropriate  values  to  use  since  there  are  no  direct 
measurements  of  the  Arctic  haze  giant  particles.  If  the  modeling  efforts 
do  show  that  the  infrared  effects  are  significant,  then  better  determina¬ 
tions  of  the  infrared  properties,  absorption  properties  in  particular, 
will  be  needed. 
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5.  CONCLUSIONS 


We  have  analyzed  the  polar  nepheloaeter  data  to  determine  aerosol  scat¬ 
tering  cross  sections.  We  have  also  used  scattering  and  aerosol  microphy- 
sical  data  to  Improve  estimates  of  aerosol  properties  for  radiative  model¬ 
ing.  We  have  confirmed  earlier  estimates  of  the  submicron  size  distribu¬ 
tion  and  have  shown  that  size  measurements  are  consistent  in  indicating  a 
mass  peak  at  ^  0.3  ym  diameter.  The  size  distributions  appear  to  be 
adequately  represented  by  log-normal  distributions  with  a  standard  deviation 
of  ^  1.5.  The  data  in  general  indicate  the  importance  of  the  submicron 
accumulation  mode  in  producing  optical  effects  at  solar  wavelengths. 

We  have  shown  that  the  giant  particles  as  measured  during  MONEX  have 
relatively  uniform  distributions  characterized  by  a  volume  peak  at 
^  2.5  ym  diameter,  a  size  somewhat  smaller  than  previously  inferred  by 
Patterson  et  al. 

We  have  also  used  a  comparison  of'  several  data  sets  to  indicate  some 
inconsistencies  in  the  data  and  to  discuss  the  reasons  for  the  inconsis¬ 
tencies.  We  have  been  able  to  point  out  some  deficiencies  in  the  optical 
particle  counter  data. 

We  have  plans  for  future  work  with  the  nephelometer  data  to  explore 
AGASP  data  in  more  detail.  For  the  purposes  of  supporting  the  conclu¬ 
sions  of  this  report,  however,  we  wish  only  to  point  out  our  new  approach 
for  analyzing  the  polar  nepheloemter  data  appears  to  offer  a  considerable 
amount  of  promise  for  future  aerosol  studies.  It  clearly  offers  a  more 
accurate  procedure  for  establishing  an  absolute  value  for  the  optical 
cross  section  of  the  aerosol  particles  by  eliminating  problems  associated 
with  trying  to  establish  the  local  number  density  to  remove  the  Rayleigh 
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scattering  cross  section  from  each  observation  of  the  atmospheric  volume 
scattering  cross  section.  Furthermore,  the  ability  to  establish  an  accurate 
value  of  the  aerosol  scattering  contribution  allows  one  to  calculate  the 
ratio  of  aerosol  to  molecular  scattering  which  was  one  of  the  variables 
obtained  in  the  least-square  curve  fitting  procedure  described  by  Grams 
(1981).  The  elimination  of  the  need  to  produce  an  estimate  of  this 
variable  will  significantly  increase  the  speed  and  reliability  of  the  pro¬ 
cedure  that  we  used  in  the  past  for  analyzing  the  nepheloueter  results. 
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